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Detecting dark matter as it streams through detectors on Earth reUes on knowledge of its phase 
space density on a scale comparable to the size of our solar system. Numerical simulations predict 
that our Galactic halo contains an enormous hierarchy of substructures, streams and caustics, the 
remnants of the merging hierarchy that began with tiny Earth mass microhalos. If these bound or 
coherent structures persist until the present time, they could dramatically alter signatures for the 
detection of weakly interacting elementary particle dark matter (WIMP). Using numerical simula- 
tions that follow the coarse grained tidal disruption within the Galactic potential and fine grained 
heating from stellar encounters, we find that microhalos, streams and caustics have a negligible 
likelihood of impacting direct detection signatures implying that dark matter constraints derived 
using simple smooth halo models are relatively robust. We also find that many dense central cusps 
survive, yielding a small enhancement in the signal for indirect detection experiments. 

PACS numbers: 98.80 



Introduction 

In a ACDM dominated universe, structure forms by the 
hierarchical clustering and merging of small density per- 
turbations [T^. Numerical simulations that follow these 
processes predict that our Galactic halo should contain a 
vast hierarchy of surviving substructures - the remnants 
of the entire halo merger tree [2 . The number density 
of substructures of a given mass M goes as n ex M~^^^ 
and they span over 15 decades in mass [3]. The small- 
est, oldest and most abundant are Earth-mass microha- 
los with a half mass radius of 10~^ pc that formed at 
2; 80 — 30 ^4 - 8 . This minimum mass is modulated by 
the free streaming velocity which is related to the mass 
of the neutralino [H [10] . 

Simulations of relatively large subhalos suggest that 
their gravitational interactions with a disk potential, can 
lead to a destruction of subhalos at distances closer than 
30 kpc [11 . Smaller subhalos form earlier, however, with 
denser cores, and are therefore the most probable dark 
matter structures to survive gravitational interactions. 
A second source of fine grained structure to survive are 
the numerous caustic sheets and folds that form due to 
the very high initial phase space density of the cold dark 
matter particles [12j. These are wrapped in a complex 
way within all the subsequent structures that form, how- 
ever in the absence of a heating term, the fine grained 
phase density would be preserved. 

With low internal velocity dispersion and high mean 
density, both the event rate and the characteristic spec- 
trum of energy deposited by dark matter in direct de- 
tection experiments [13] could be affected by any fea- 
tures surviving in the phase space distribution of CDM 
particles. Direct detection experiments are sensitive to 
the density and velocity distributions of WIMPs on a 
scale of ^ lO^^m, the distance the Earth travels over a 
year. In order to make predictions and exclusion lim- 



its, these experiments assume that the dark matter is 
completely smooth on these scales, with a well mixed 
Maxwell-Boltzmann velocity distribution [Ml |T5] . Fur- 
thermore, if any such small high-density clumps actually 
dominate the dark matter distribution in the solar neigh- 
borhood, the indirect detection signal due to dark matter 
annihilation in the galaxy might also be affected. 

Since existing N-body simulations of galaxy formation 
do not have a resolution that goes down to objects with 
mass as small as IO^^Mq, in order to address the ques- 
tion of the survival and impact of such microhalos we 
need to combine analytical estimates with the results of 
smaller scale simulations that can resolve such objects. 

Previous studies have made analytic [16 and numerical 
estimates [17] of the disruption timescale of microhaloes 
as they orbit through the stellar field. Zhao et. al. [16] 
argued that most of the microhaloes should be completely 
destroyed by encounters with stars, whilst Goerdt et. al. 
[17] show that indeed, whilst most of the mass is unbound 
the dense central cusp may survive intact. 

We extend this work by numerical calculating the tidal 
disruption of microhaloes as they actually orbit through 
a field of stars, as well as self-consistently including the 
Galactic halo potential. We calculate the survival statis- 
tics of microhaloes using realistic orbital distributions 
within the disk, allowing us to follow the dynamical struc- 
ture of the dark matter streams and thus to estimate the 
fine grained phase space distribution function of WIMPs 
on scales relevant to dark matter detection experiments. 



Microhalo parameters and disruption processes 

The dominant processes that can affect microhalos in- 
volve gravitational interactions with baryons in the stel- 
lar field during the crossing of the disk and also tidal 
effects of the disk potential during the orbit of the micro- 
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halo. Unfortunately we cannot account for both effects 
simultaneously in our simulations as it would require set- 
ting up a self-consistent disk with billions of stars. In- 
stead we look separately at the effects of stellar disruption 
and tidal streaming and we then estimate the combined 
behavior. 

As a substructure halo crosses through a stellar field, 
high-speed interactions with single stars will heat up the 
halo distribution, causing it to increase its velocity dis- 
persion and hence its scale size will grow. This process 
is analogous to galaxy harassment that occurs in clus- 
ters [T8^ and basically has a timescale proportional to 
the relaxation time of the stellar disk and the time each 
microhalo spends within the fluctuating potential field of 
the stars. For an analytical estimation of this process we 
follow the Goerdt et. al. paper [17]. 

In the 'distant-tide' approximation [19] the internal en- 
ergy increase of the microhalo due to a single encounter 
with a fixed star is given by 
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where b is the impact parameter and Vmh is the velocity 
of the microhalo. Since halos with an early formation 
time have a low concentration we can set (r) ~ O.Sr^^^. 

One encounter can totally disrupt the microhalo if 6E 
exceeds the binding energy E^j. Since E^j ^ OAv^^^ with 
'^vir — Grriyir/ryir [20 , the minimal encounter parameter 
that does not entirely disrupt the microhalo is found to 
be 
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We can now define the disruption probability of a micro- 
halo in a stellar field 



p = ^ I SEdN 



dN^— SEdN, (3) 

Jbrr^^r^ 



where dN = 2t: nhdhVmhdt. Here we have used 5E = 
El) for h < bmin- Performing the integration leads to a 
disruption probability of 
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where we have used the definition of the virial radius 
Myir = ^r^^^200pc with pc = The microhalo is 

completely destroyed at p = 1. Therefore we get the 
average disruption time 



t = 250 
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Myr. (5) 



A microhalo with a formation redshift z 60 should 
therefore survive about 250 Myr in a stellar field with 
a density similar to the one in the solar neighbourhod. 



However this is only true on average, since one very close 
encounter can immediately lead to total disruption. 

The above estimate does not take into account the in- 
ternal structure of the microhalos and should therefore 
only give a very rough estimation of the disruption time. 
Also it does not allow one to follow the mass loss during 
the disruption process. We therefore perform a simula- 
tion where a microhalo is crossing a periodic box of stars. 
The box has a length of 50 pc and is filled with randomly 
distributed stars with the density p = O.O4M0pc~^ and 
the velocity dispersion a = 50 km/s. This constellation 
corresponds to the stellar field in the disk at the solar 
radius [19 . For simplicity all the stars have the average 
mass of 0.7 M©. The microhalo which is crossing the box 
at 200 km/s is set up by the halogen-code of Zemp et. 
al. [21 . Corresponding to the results in [7] it has a mass 
of 1O~^M0 and the density profile 
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with a = 1, /3 = 3 and 7 = 1.2, as well as a concentration 
parameter of c = r2oo/^s = 1-6. The virial radius r2oo is 
defined with respect to the background density at 2: = 60, 
the average formation redshift of a microhalo [4 . 

In our simulation we find that 50% of the microhalo 
mass is unbound after 80 Myr of box-crossing (which 
corresponds to about 40 perpendicular disk passages). 
After 160 Myr (80 disk passages) even the central core 
starts to disappear and more than 90% of the microhalo 
is completely disrupted (see pictures in Table [ij). At latest 
after 200 Myr (100 disk passsages) no bound structure is 
left (see Fig[T]). Our simulation gives therefore a slightly 
shorter disruption time than the simplified analytical es- 
timation of equation ([5|. 

In order to determine what fraction of microhalos sur- 
vive until the present day, we have to calculate orbital 
statistics and the distribution of disk crossing times. This 
can be established by tracing back the orbits of particles 
in a galactic potential. We use the standard Milky Way 
model with disk and halo particles set up by the Galac- 
tlCS code [22] and we select a sample of halo particles in 
a small box around the position of the sun. The orbits 
of these particles are followed backwards in time and we 
find that the average number of disk crossings for these 
particles is c = 80 with a standard deviation of dc = 43. 
The average crossing radius is (not surprisingly) R = S 
kpc with (jR = A kpc. The spread of disk crossing events 
for different particles follows a Maxwell-Boltzmann dis- 
tribution. 

We use this disk crossing distribution combined with 
the rate of mass loss determined from our numerical 
study to calculate the survival statistics of microhalos 
in the vicinity of the sun. Since the timescale for com- 
plete disruption in our simulation is equivalent to the 
average time a microhalo spends in the stellar disk, we 
conclude that the average microhalo in the vicinity of the 
sun is just about to be entirely destroyed at the present 




TABLE I: Microhalo density map att= 0, 20, 40, 60, 80, 100, 
120, 140, 160 Myr (from the upper left to the lower right). 
The boxlength of the images is 0.38 pc. 

time (see also [23]). At most five percent of its initial 
mass is still in a bound core. However the spread in the 
number of disc crossings is relatively wide and a signif- 
icant fraction of microhalos should still have surviving 
cores. Mass loss is nevertheless important: microhalos 
maintaining more than 50% of their initial mass should 
be rare. Figure [l] illustrates the mass loss, where the red 
curve shows the disruption of a typical microhalo with 
80 disk crossings in 10 Gyr at the radius of the sun. 

However, disk crossing is not the only source of dynam- 
ical disruption. While orbiting the galaxy, a microhalo 
is under the constant influence of the global Galactic po- 
tential, and tidal forces will act so that the microhalo's 
structure becomes elongated and unbound particles will 
form leading and trailing tidal streams. The detailed 
impact of tidal streaming depends on the orbit of the mi- 
crohalo and on the shape of the host potential. In our 
simulations we use a disk potential that emerges from a 
density distribution of the form 

p{R, z) (X exp{-R/Rd)sech^{z/zd). (7) 

Here R and z are the disk radius and the height respec- 
tively, which we set to be Rd = 2.8 kpc, Zd = 0.4 kpc. 
The disk mass is = 4.5 • 10^^ Mq. In all our simula- 
tions the orbit of the microhalo is chosen to be roughly 
spherical with a distance of 7.9 kpc from the galactic 
center. 

We cannot model both heating due to stellar inter- 
actions and tidal elongation at the same time since this 
would require following the motion of 50 billion disk stars. 
We therefore performed orbital simulations for three dif- 



FIG. 1: Ratio between bound and total mass of a microhalo 
crossing a stellar field (red) and orbiting in a Milky Way po- 
tential (black, grey). The red curve is stretched in order to 
simulate the effect of disk crossing on an average microhalo 
with about 80 disk crossings. The black and the grey curve cor- 
respond to a microhalo that spent Myr respectively 80 Myr in 
the stellar field before orbiting around the galactic potential. 
The dominant mass loss is coming from stellar interaction 
until complete disruption after about 12 Gyr. 



ferent cases: an initially completely undisturbed micro- 
halo, a microhalo that first crossed the stellar field for 
80 Myr and has lost about 60 percent of its mass, and 
a completely disrupted microhalo that spent more than 
160 Myr in the stehar field. 

The length of the tidal streams / due to the orbit- 
ing process can be crudely estimated with the relation 
l{t) ~ cFmh^i where (Jmh is the velocity dispersion of the 
initial microhalo. For the initially unperturbed microhalo 
(Jmh ^ 10~^ km/s, causing a stream length of roughly 
/ ~ 10 pc after one Hubble time. For the initially com- 
pletely disrupted microhalo Gmh ^ 10~^ km/s, and the 
stream length is about / ~ 100 pc after a Hubble time. 
These length scales agree well with our simulations (i.e. 
see Table [ll|. 

Orbiting in the galactic potential significantly reduces 
the mass of the microhalo (see black and grey lines in Fig 
[1]). However, the rate of tidal mass loss is suppressed as 
the tidal radius is steadily reduced. The central cusp of 
each dark matter microhalo has a very deep potential, as 
a consequence there is always a bound core remaining, 
even for a microhalo that has been heated in the stellar 
field before orbiting. 

Comparing the curves in Fig [l] leads to the conclusion 
that disk crossing is the dominant disruption process and 
the only one that can lead to complete distruction of 
the microhalo. The step-like decrease of the curve is an 
indication of very close encounters that play a mayor role 
in the disruption process. Tidal stripping on the other 
hand can also significantly reduce the mass but it never 
completely distroys the microhalo because of its tightly 
bound inner core. 



TABLE II: Streaming microhalo after 10 Gyr on a roughly 
circular orbit around a Milky Way potential: The two images 
on the top show the sheet-like streams from the top and from 
the side (boxlength: 30 pc). The third image is a zoom in at 
the centre where the still hound core is visible (boxlength: 2 
pc). For these pictures we have used an initially non disrupted 
microhalo with formation redshift 60. 



Implications for Dark Matter Detection 

In direct detection experiments the differential inter- 
action rate is sensitive to the fine grained density and 
the velocity distribution of dark matter particles on A.U. 
scales [241 123 • Substructures like microhalos can affect 
the interaction rate if they are abundant enough to have 
a substantial likelihood of existing in the solar neigh- 
bourhood and if their density is at least the same order 
of magnitude as the background dark matter density in 
this region, pi^g ~ lO^M0kpc~^ (see for example 
Equally important, the phase space for the energy de- 
posits associated with dark matter events will not be that 
appropriate for an isothermal halo if a single microhalo 
were to dominate the density distribution in the solar 
neighbourhood [T3[ [27] . 

Our results above suggest that none of these condi- 
tions are generally achieved. In Figure [2] we plot the 
stream densities of microhalos that crossed the stellar 
field for Myr (red), 80 Myr (grey) and 160 Myr (blue), 
before orbiting in the galactic potential for 10 Gyr. The 
tidal streams of the initially unperturbed halo (red) have 
an average density oi p ^ lO^M0kpc~^, which is al- 
ready negligibly low compared to the background. Only 
the very tiny core still maintains its initial density of 
p ^ lO^^M0kpc~^. The initially disrupted microhalo 
(blue) has no more bound core. Its stream density is 
only at about p ~ 10^ — lO^M0kpc~^. The approach 
of first measuring the stellar disruption and then looking 
at tidal effects underestimates the stream density, since 
the microhalos get most of their heat energy right at the 
beginning. Therefore, the actual stream density of an 
average microhalo should be somewhere between the red 
and the blue line in Figure [2] 

Since the stream densities are far below the value of 
the local galactic density, only a surviving core existing 
in the region of the earth would any effect upon direct 



r [kpc] 

FIG. 2: Stream densities of microhalos after an orbital time of 
ten Gyr. Before orbiting the microhalos have spent a time of 
Myr (red), 80 Myr (grey) and 160 Myr (blue) in the stellar 
field. There is still a visible bound core in the red and the grey 
profile but no more in the blue one. The black dots represent 
the density profile of a completely undisrupted microhalo that 
spend no time on an orbit. The straight blue line corresponds 
to the average dark matter density around the sun. 



detection. However, only about half of the microhalos 
still have bound cores because of disk crossing, and tidal 
effects further reduce the mass of the cores to less than 
ten percent of their original value. We note that any sub- 
structures orbiting primarily within the disk plane would 
be quickly destroyed by stellar encounters. 

The chance of being in such an overdense region can 
then be optimistically estimated: An extrapolation of 
the subhalo mass function leads to a microhalo number 
density rimh of about 500 pc~^ at the solar radius 0. 
This number can be divided by two due to the disruption 
processes stated above and again by two since microhalos 
orbiting in the disk plane are completely disrupted. We 
then end up with the approximation of rimh ^ 100 pc~^ 
at the solar radius. Each microhalo has a volume of about 
Vmh ^ 10~^ pc^ and therefore there is a chance of about 
0.0001% of being in such an overdense region. 

The streams of particles stripped from microhalos are 
coherent and long, thus it is appropriate to calculate 
their volume filling factor. Since the stream density is 
p ~ 10^ — lO^M0kpc~^ we expect that our solar sys- 
tem is criss-crossed with /5 x (10^ — 10^) streams, where 
fi) ~ 0.1 is the fraction of the local Galactic halo den- 
sity that forms from substructures up to a solar mass. 
Larger substructures may be completely disrupted at the 
Sun's position in the Galaxy due to global disk shocking 
and tides [11] . The velocity dispersion within an average 
stream due to heating by disk stars is a ~ 10~^ km/s. 
Thus, the local density is determined by the superposi- 
tion of a large number of independent streams, and the 
overall velocity distribution at the solar radius should 
be essentially Maxwellian, isotropic and smooth with no 
spiky structure, as we would assume for a smooth halo 
model with no substructures. The signatures of streams 
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could be only be detected experimentally with over sev- 
eral hundred events. 

The case for indirect detection is somewhat differ- 
ent from that described above. In indirect detection 
experiments one tries to detect the annihilation prod- 
ucts, such as gamma-rays, coming from the highest den- 
sity dark matter regions, which is proportional to the 
square of the dark matter density times the volume of 
the region observed ^28H31j . Consider a volume con- 
taining on average one microhalo V ~ 10~^ pc~^. The 
luminosity due to the smooth background is therefore 
Lbg oc Vp'^ = 10~^MqPc~^, whereas the luminosity of a 
surviving microhalo core is 



Lmh OC Vcoreplore ^ 5 X 10 ^ M^pC 



^-3 



(8) 



where we have assumed a mean core density of 
IO^^MqPc"^. Thus the net boost factor due to micro- 
halos is about 1.5, and stays below the detection limits 
of the FERMI experiment [32]. However this number 
is highly uncertain since it depends on extrapolations of 
both the substructure mass function and also on the mi- 
crohalo internal density structure. 

Finally, we consider the evolution of caustic sheets of 
particles within the Galactic halo. In the absence of fine 
grained heating, narrow sheets and folds will occupy re- 
gions of phase space within all collapsing CDM struc- 
tures. As structures merge hierarchically, these caustic 
features become wrapped in phase space like a fine fab- 
ric that has been crumpled into a ball, the phase space 
density at any point being preserved. During the mat- 
ter dominated epoch and before structure formation the 
velocity dispersion of WIMPs is given by 



_ in flOOGeVV^^ 
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where is the mass of the WIMP 



. Since the first 
structures are collapsing at redshifts z ~ 60, we obtain 
a primordial velocity dispersion of ~ 2 cm/s. In the 
outer halo these features will persist, but in the vicinity of 
the sun heating by the disk stars will broaden the phase 
space and physical space distribution. 

Since the energy increase is proportional to time (see 
equation ([3|) we obtain a heating factor 
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where C has been determined via our simulations. For 
a microhalo on an average orbit we find an increase in 
the velocity dispersion of Aa ~ 12 m/s which is much 
larger than the primordial velocity dispersion, effectively 



smearing out all caustic overdensities in the vicinity of 
the Galactic disk. 



Conclusion 

As the prospects for direct and indirect detection of 
WIMP dark matter improve with the development of 
new detectors, a renewed interest in the phase space 
distribution of dark matter particles has arisen. It has 
recently been shown that hierarchical clustering contin- 
ues down to extremely small mass scales, so that most 
dark matter currently in the halo of our galaxy may 
have originated in microhalos with masses as small as 
1O~^M0. Possibly dense surviving cores, tidal streams, 
and caustic structures might leave phase space sparsely 
populated, suggesting exciting new possibilities for novel 
signatures that differ from the traditional experimental 
assumption a smooth isothermal halo. However, our re- 
sults imply that tidal effects and gravitational heating 
effectively wipe out any such signatures for Earth based 
detectors. Even though we find that a significant fraction 
of microhalos still have a bound core today, these over- 
dense regions are too small to be relevant for detection 
experiments. 

The disrupted material in the tidal streams is not dense 
enough to affect the detection signal. On Earth, there are 
about 7 • 10^ dark matter particles per second streaming 
through our bodies (assuming ~ 100 GeV), but they 
originate from over 10^ streams coming from disrupted 
microhalos. The velocity dispersion in the streams is 
heated up through stellar interaction from initially 10~^ 
km/s to 10~^ km/s and therefore we expect an essentially 
smooth Maxwellian, with at most some spikes due to 
microhalos with an unusual orbital history [33l |34] . 

To summarise, our results imply that limits obtained 
on dark matter from detection experiments under the 
conservative assumption of a smooth halo with nearly 
Maxwellian density distribution remain valid, with a 
prefactor depending only on the average local dark 
matter density. The characteristic deviations from a 
Maxwellian distribution predicted from numerical simu- 
lations may be detected given sufficient detection statis- 
tics [35 , 36 . Substructure cores still persist in the vicinity 
of the Sun to today, however the boost factor for indirect 
detection even for these is likely to be relatively small. 
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